Voter Models on Heterogeneous Networks 
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We study simple interacting particle systems on heterogeneous networks, including the voter 
model and the invasion process. These are both two-state models in which in an update event an 
individual changes state to agree with a neighbor. For the voter model, an individual "imports" its 
state from a randomly-chosen neighbor. Here the average time Tm to reach consensus for a network 
of N nodes with an uncorrelated degree distribution scales as Nfii/fi2, where fxj, is the k th moment 
of the degree distribution. Quick consensus thus arises on networks with broad degree distributions. 
We also identify the conservation law that characterizes the route by which consensus is reached. 
Parallel results are derived for the invasion process, in which the state of an agent is "exported" to 
a random neighbor. We further generalize to biased dynamics in which one state is favored. The 
probability for a single fitter mutant located at a node of degree k to overspread the population — the 
fixation probability — is proportional to k for the voter model and to 1/k for the invasion process. 
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I. INTRODUCTION 

Recent studies of statistical physics models on com- 
plex networks have elucidated the effect of heterogene- 
ity in link structure on dynamical properties and critical 
behavior. For scale-free networks, the source of hetero- 
geneity is the broad distribution of node degrees, where 
node degree is defined as the number of links attached to 
a node. This dispersityleads, for example, to a vanishing 
percolation threshold [1| , allows epidemics to thrive even 
with a vanishingly small infection rate [2j , and causes the 
Ising model to be ordered at all temperatures 0]. 

For many of these models, the dynamics can be un- 
derstood by accounting for the broadness of the degree 
distribution within a mean-field description. In this arti- 
cle we show how to implement such an approach for two 
of the simplest interacting particle systems on heteroge- 
neous networks, namely, the voter model (VM) 0, Q and 
its close relative the invasion process (IP) [6| , as well as 
their generalizations to biased evolution |U [l(| • 

For the VM and the IP, each node of an iV-node net- 
work can be in one of two discrete states: and 1. In a 
social context, these states represent two possible opin- 
ions of the individual at that node. In a biological con- 
text, the states represent the phenotype of the individual, 
with state representing resident individuals and 1 for 
mutant individuals. The evolution consists of changing 
the state of a node at a rate that depends on its local 
environment. The difference between the VM and IP is 
in the order in which the "invader" node (whose state is 
adopted) and the "receiver" node (whose state changes) 
are chosen. This ordering is immaterial on degree-regular 
graphs, such as regular lattices, but it plays an essen- 
tial role on heterogeneous degree networks. Understand- 
ing the basic difference between the VM and IP on such 
graphs is one of the main goals of this work. 

We will also study the biased voter model and the bi- 
ased invasion process in which there is a preference for 



one of the two states. This generalization describes the 
evolution of a fitter mutant in an otherwise homogeneous 
resident population. We will determine the probability 
for a single fitter mutant to overspread a population with 
biased VM and biased IP dynamics. Our main result 
here is that the probability for a single fitter mutant at 
a node of degree k to overspread a population — the fixa- 
tion probability — is proportional to k for the VM and to 
1/k for the IP. 

In Sec. HH we define the models of this paper. In 
Sec. IIII1 we summarize the properties of the VM on the 
complete graph. We then investigate the VM on het- 
erogeneous networks in Sec. IIV1 including the complete 
bipartite graph, the two-clique graph, and general degree- 
heterogeneous graphs. In Sec. [Vj we treat the comple- 
mentary IP. Finally, in Sec. lVIl we investigate the biased 
versions of the VM and the IP and determine the fixation 
probabilities on general heterogeneous networks. 



II. FORMULATION OF THE MODELS 



A. Update Rules 



VM evolution consists of the following two steps: 

(i) pick a random node (a voter); 

(ii) the voter adopts the state of a random neighbor. 

In the closely related IP the evolution steps are: 

(i) pick a random node (an invader); 

(ii) the invader exports its state to a random neighbor. 

We also mention an intermediate model, link dynamics 
(LD), whose evolution steps are: 

(i) pick a random link; 

(ii) one of the nodes on the link, adopts the state of the 
other end node. 
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In these models, the time is incremented by 1/JV in each 
update. Thus JV updates corresponds to each node being 
updated once, on average, after which the time increases 
by 1. Steps (i) and (ii) are then repeated ad infinitum 
or until the system reaches consensus, an event that is 
certain to occur in a finite time when the network is finite. 

pick 1st pick 2nd 



(a) 



pick 1st pick 2nd 



(b) 



=0^ 1- 



pick link 



=(k — 1 




(c) 



FIG. 1: Illustration of the update rules, (a) Voter model 
(VM): a voter is chosen at random and it adopts the state 
of random neighbor, (b) Invasion process (IP): a randomly 
chosen voter exports its state to a neighbor, (c) Link dynam- 
ics (LD): an active link is randomly picked and a randomly 
chosen node at the end of this link is updated. 

The interpretation of VM dynamics is that individu- 
als lack any self-confidence. In an update, an individual 
therefore consults one of its neighbors and adopts the 
neighbor's state. On the contrary, in the IP an individ- 
ual imposes its state to one of its neighbors. Here we can 
think of the selected individual as replicating, and its 
offspring invades and replaces the individual at a neigh- 
boring node. While the different update details of the 
three models might appear superficially trivial, we shall 
show that these differences are fundamental when the 
dynamics occurs on degree-heterogeneous networks. 

In general, such a network may be specified by its ad- 
jacency matrix A, with A xy = 1 if nodes x and y are 
connected, and A xy = otherwise. The degree distribu- 
tion of such a heterogeneous network is specified by 



n k 



Nk 
N 



(1) 



where JV& is the number of nodes of degree k and JV is 
the total number of nodes in the network. The moments 
of the degree distribution are then 



n k . 



(2) 



Special cases of relevance for the VM and IP are the av- 
erage degree /ii, the second moment fa, and the average 
inverse degree fi—i. Normalization also fixes the zeroth 
moment /io = 1. 

Let 77 represent the state of the entire network, and 
define r){x), which can take the values or 1, as the state 
of node x. In each update event in the three models, the 
state of a single node changes from to 1 or vice versa 
(Fig. [T]) . We represent by r\ x the state of the system that 
results after changing the state of the node at x: 



Vx(y) 



r](y) y^x 
1 — f](x) y = x 



(3) 



Transitions are specified by the probability that the state 
of node x changes, which we term a flip event. For the 
VM, IP, and LD, the transition probability at node x is 



.4 



(4) 
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LD, 
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where $(x,y) = r\(x)\Y — T](y)], and the quantity 

A xy [<f>(x iy ) + <P{y,x)} 

is non-zero only when nodes x and y are connected (the 
factor A xy ) and in opposite states (expression in square 
brackets), so that an update can actually occur. Finally 
Q is 



(5) 



For the VM, the factor (JVQ) -1 in (gj) accounts for first 
choosing node x with probability 1/JV, and then any of 
its neighbors y with probability 1/Q = l/k x . Conversely, 
in the IP, one first chooses node y (a neighbor of x) with 
probability 1/JV, and then one chooses x with probability 
1/Q = 1/ky. In LD, each link is chosen with the same 
probability 2/JV/xi and then one of the nodes at the end 
of this link is updated with probability 1/2, which then 
leads to Q — /xi. As we shall discuss, VM, IP, and LD 
dynamics are distinct for networks with heterogeneous 
node degrees. 

The above models have been defined as discrete-time 
processes, where one node is chosen in each update 
step. Alternatively, these models can be formulated as 
continuous-time processes by increasing the time by a 
random increment that is chosen from an exponential 
distribution with mean 1/JV. This stochastic update is 
equivalent to each update event occurring with unit rate 
in the three models. For both discrete and continuous 
time dynamics, the fixation probabilities and the average 
time to consensus (for long times) are the same. How- 
ever, there is an essential difference between discrete and 
continuous time dynamics for biased models that will be 
discussed in Sec. IVTl 
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B. Conservation Laws and Exit Probability 

An important aspect of the VM, IP, and LD is that 
each model has its own dynamically conserved quantity. 
This conservation law determines the fundamental exit 
probability £i(p), namely, the probability that a finite 
system with an initial density p of Is reaches a consen- 
sus of all Is. This quantity is also called the fixation 
probability in the biology literature. 

The simplest case is LD, for which the state of the sys- 
tem changes only when an active link — where the nodes 
at the ends of the link are in opposite states — is cho- 
sen. Because the invader and the receiver are assigned 
randomly, the probability of increasing or decreasing the 
number of 1 nodes (mutants) are the same. The dynam- 
ics is thus equivalent to a symmetric random walk on the 
integers, with absorbing states at and at N mutants. 
Because of the symmetry of the random walk, the ini- 
tial and final densities of nodes in state 1 are the same. 
Consequently, the exit probability is 



where 



£i(p) = P = Jj^2v(x)- 



(6) 



This result is also valid for the VM and the IP on degree- 
regular graphs, due to their equivalence to LD when every 
node has the same degree. 

We now extend Eq. © to general networks. Consider 
the average change in n{x) at node x, (Ar](x)}. Here the 
angle brackets denote the average over all realizations of 
the update dynamics. This change equals the probability 
that rj(x) increases from to 1 minus the probability that 
r](x) decreases from 1 to 0. Hence 



(A v (x)) = [l-2»j(x)]P [v^ifc] 



(7) 



Substituting the transition probability from Eq. ([J), we 
obtain 



(8) 



where we use the fact that rj(x) 2 = rj(x). The change 
in the average density in the entire network, (Ap), is 
obtained by summing over all nodes x to give 

(Ap) = "£{A V (x)) = £ || foGO ~ »?(*)] ■ ( 9 ) 

x x,y 

For LD, and for any of the three models on regular 
graphs, Q is constant. Hence the summand in the ex- 
pression on the right is antisymmetric in x and y and 
(Ap) =0. As a consequence (p) is conserved, so that 
the mutant fixation probability equals the initial mutant 
density p. 

To compute the exit probability on arbitrary networks 
for general models, we generalize the notion of density 
by introducing the degree-weighted moments 



1 



Np r , 



£«*) 



i 

f-^rn 



k 



n k pk 



(10) 



Pk 



(11) 



is the density of Is on the subset of nodes of degree k. 
Here the prime on the sum denotes the restriction that all 
nodes x have fixed degree k. Note that ujq coincides with 
the density p. To obtain a conserved quantity, it is clear 
that the factor Q in the denominator of the transition 
rate in must be canceled out. For the VM, we thus 
consider A(a>i) and repeat the calculation that led to 
Eq. © to obtain 



A ^>=E]v^ [r){y)-r>(x)\ 



Similarly, the IP, we consider A(u>_i) and obtain 



A( W _ 1 ) = E 



-4, 



(12) 



(13) 



Because of the x-y antisymmetry of the summand in 
Eqs. (fT2"]) and (|13p. both sums vanish so that the con- 
served quantity in the three models are: 



("i> 

(w ) = (p) 
(w-i) 



VM, 

LD, 

IP. 



(14) 



We can understand these conservation laws intuitively. 
For the VM, although nodes are selected uniformly, there 
are relatively more low-degree nodes that are neighbors 
of high-degree nodes (see Fig.[T]). Thus low-degree nodes 
change their state more often than high-degree nodes. 
Weighting each node by its degree compensates this dis- 
parity and leads to the conserved quantity u>i . Thus the 
mean density p is not conserved, as first pointed out in 
Ref. [ll| for the VM on heterogeneous graphs. In con- 
trast, For the IP low-degree nodes change their state less 
often than high degree nodes, and this disparity may be 
compensated by weighting each node by its inverse de- 
gree. This leads to the conserved quantity uj-\. 




FIG. 2: A star graph. 

Since the initial value of the conserved quantity equals 
its value in the final unanimous state, we obtain, for the 
exit probability 



£i(<*>o) =ujq = p 



VM, 

LD & degree regular graphs, 
IP, 



(15) 
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An instructive example is the extreme case of a star 
graph, where N nodes are connected only to a single 
central hub (Fig. [2]) . For the VM, if the hub is in state 1 
and all other nodes are in state 0, then Eq. (fT5"|) predicts 
that the probability of reaching 1 consensus is 1/2! Thus 
a single individual with a macroscopic number of neigh- 
bors largely determines the final state. Conversely, for 
the IP, this same initial state reaches 1 consensus with 
probability that is 0(1/N 2 ). We will discuss this dra- 
matic disparity between VM and IP dynamics in detail 
in Sec. ED 



III. VOTER MODEL ON THE COMPLETE 
GRAPH 

As a preliminary for degree-heterogeneous graphs, we 
discuss the well-known dynamics of the VM on the com- 
plete graph 0, [12, HH , where each node is connected to 
every other node. Because the complete graph is degree 
regular, the VM, IP, and LD are all equivalent and we 
treat the system in the framework of the VM. Let p(t) 
be the density of voters in state 1. In each update event 
p —>■ p±5p, with Sp — 1/N, corresponding to the respec- 
tive state changes — > 1 or 1 — > 0. The probabilities for 
these events are 



R(p) EE P[ P 



■Sp] = (l-p)p 



L(p) = P[p^p-Sp]=p(l-p). 



(16) 



Here R and L denote raising and lowering operators that 
give the transition probabilities from p to p± dp, respec- 
tively. 

Let c(p,t) be the probability that the density of l's is 
p at time t. After one update event, this density evolves 
according to 

c(p, t + 8t) = R(p — Sp)c(p — Sp,t) 
+L(p + 5p)c(p + 8p,t) 
+[1-R(p)-L(p)]c(p,i). (17) 

Here St = 1/N and the first two terms on the right ac- 
count for the inflow to the state with density p in an 
update and the last term accounts for outflow. Expand- 
ing Eq. ([17)) to second order in Sp gives the forward Kol- 
mogorov or Fokker-Planck equation [l4| , 



dc(p,t) d 



d 2 



uj =—[ v ( p )c(p,t)} + —[D(p)c(p,t)}, (18) 



where 



Sp, 



v{p) = ^mp)--L{p)\, 

is the drift term (called the selection term in biology (l5j ) 
that is caused by the bias in the transition probabilities, 
and 

D( P )e1|[R(p) + LW], 



is the diffusion term (paradoxically called the random 
drift term in biology) that quantifies the stochastic noise 
in the kinetics. On the complete graph, the selection term 
is zero and the Fokker-Planck equation ([T8"|) becomes 



dc{p,t) _ 1 3 



at 



Nd fW-PHP,t)]. 



(19) 



In a similar fashion, the equation for the exit proba- 
bility with initial density p is 

£ 1 (p) = R(p)£ 1 (p+8p) + L(p)S 1 (p-8p) 

+ [1-R(p)-L(p)]£ 1 (p). (20) 

This equation expresses £\ (p) as the probability of mak- 
ing a transition to p±8p or p, respectively, times the exit 
probability from this intermediate point [131 ] . Expanding 
Eq. ([2U)l to second order in Sp gives the backward Kol- 
mogorov equation for the exit probability 



v(p) 



d£i(p) 
dp 



D{ P ) 



d 2 £!( P ) 
dp 2 



B£i = 0, (21) 



with B the generator of the backward Kolmogorov equa- 
tion. For the boundary conditions £\{p = 0) = and 
£i{p — 1) = E the solution is simply £\{p) — p. This 
result reproduces Eq. (jl5| that was obtained by magne- 
tization conservation. 

In analogy with Eq. (|20|) . the average time to reach 
consensus, T(p), as a function of the initial density p, 
obeys the backward Kolmogorov equation [131 ] 

T(p) = St + R(p)T(p+Sp) + L(p)T(p-S P ) 

+ [1 - R(p) - L(p)]T(p). (22) 

This equation expresses the average consensus time as 
the time for a single step plus the average time to reach 
consensus after taking this step. The three terms account 
for the transitions p — ► piSp or p — > p, respectively, and 
the factor St in each term accounts for the time elapsed 
in a single update. Expanding Eq. ([2"2"]l to second order 
in Sp gives the backward Kolmogorov equation for the 
average consensus time 



v(p)— — +D{p)- 



dp 



dp 2 



BT= -1. 



(23) 



Using the transition probabilities in Eqs. ([16|) and set- 
ting St = Sp = 1/N, the above equation reduces to 



p(l - p) d 2 T{p) 



N 



dp 2 



-1. 



(24) 



For the boundary conditions T(0) = T(l) = 0, the solu- 
tion is 



T(p) = N 



(1-p) In- +pln- 

1 -p p^ 



(25) 



which is symmetric about p = 1/2. Important special 
cases are T(l/2) = N In 2, corresponding to starting with 
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equal densities of voters of each opinion, and T(l/N) f» 
In N, corresponding to starting with a single mutant. 

In addition to the time to reach either type of consen- 
sus, consider the fixation times, namely, the conditional 
times to reach 1 consensus, defined as or con- 

sensus, To(p), as a function p. We obtain these times 
by extending the backward Kolmogorov approach [13] to 
account for the conditioning on type of consensus. The 
conditional fixation times satisfy 

B [e T ] = -£ B [£ a Ti] = -Si, (26) 

with £o = 1 — Si and B is the generator in Eq. (|23| . 
subject to absorbing boundaries for both SoTo and Z\T\. 
The solution to Eq. (gHJ) is 



Mp) = N- 



-ln- Ttip) 
P P 



N- 



.1 



In- 



1 



. (27) 

1-P P V P 1-P 

These fixation times then satisfy the sum rule 
£i(p)T x (p) + £ (p))2o(p) = T(p) that reflects the fact 
that the consensus time is the suitably weighted aver- 
age of the fixation times. One important limit is the 
initial state of a single 1 mutant, for which the fixation 
time to all Is is Ti(l/AT) « N. In contrast, the consen- 
sus time from this same starting state is much smaller, 
T(l/N) « InN, because the system can exit via the 
nearby boundary at p — 0. 



IV. VOTER MODEL ON COMPLEX 
NETWORKS 

A. The Complete Bipartite Graph 

To understand how degree dispersity affects VM dy- 
namics, we first study a simple degree-heterogeneous net- 
work, the complete bipartite graph K a ^, in which a + b 
nodes are partitioned into two subgraphs of size a and b 
(Fig. [3]). Each node in the a subgraph is connected only 
to all nodes in b, and vice versa. Thus a nodes all have 
degree b, while b nodes all have degree a. 



a nodes 
degree b 



b nodes 
degree a 



FIG. 3: (Color online) The complete bipartite graph K a ,b- 

Consider the voter model on this graph. Let A^b be 
the respective number of voters in state 1 on each sub- 
graph and let p a = N^/a, pb — Nb/b be the respective 
subgraph densities. In an update, these numbers change 
according to transition probabilities, 




Ra = P[Pa,Pb -* Pa'Pb] 
La = P[Pa,Pb -> Pa'Pb] 



Pb(l - Pa), 



a+ a h , (28) 

— T Pa 1 ~ Pb), 

a + b 



with = p a ± a -1 . Here R a is the probability to in- 
crease the number of Is in subgraph a by 1, for which 
we need to first choose a in subgraph a that then in- 
teracts with a 1 in subgraph b. Similarly, L a gives the 
corresponding the probability for reducing the number 
of Is in a. Analogous definitions hold for Rb and Lb by 
interchanging a <-> b. 

From these transition probabilities, the rate equations 
for average subgraph densities are 
flq. with solution 



dt 



Pb,a — Pa,b 



Pa,b(t) = 



Pa(0) + p b (0) [Pa,b(0)-p b ,a(0)]< 



(29) 



Thus the subgraph densities are driven to the common 
value (/O a (0) + pb(0))/2 in a time of the order of one. 
As a result, the density of Is in the entire graph, which 
evolves as 



dp 
~dt 



1 / dfa_ 

+ b\dt 



dpt 
1 dt 



b — a 



(Pa - Pb) 



becomes conserved in the long-time limit. Thus there is 
a two-time scale approach to consensus: at early times, 
there is a non-zero bias that quickly drives the system to 
equal subgraph densities p a — Pb] subsequently, diffusive 
fluctuations drive the eventual approach to consensus. 

We can understand this behavior in a more fundamen- 
tal way by studying the probability that the graph has 
density of Is equal to p a and pi, in each subgraph at time 
t, c(p a , pb, t). This probability density evolves as 

c(p a ,Pb,< + St) = 

Ra (P a . Pb ) c(Pa > Pb , t) + L a (p+ ,p h )c{p+,p h , t) 
+ R b (p a , P b ) c (Pa, P b ,*) + Lb (Pa, Pb ) c (Pa , Pb , *) 
+ [l - R a (/3a, Pb) - L a (p a , pb) 

- Rb(Pa,Pb) - L b (/9 a ,Pb)] c(p a ,pb,t). (30) 

Expanding to second order in a -1 and b^ 1 gives the 
Fokker-Planck equation, 



dc _ Id 
dt aSt dp a 
1 d 2 



2a 2 5t dpi 



[(R a -L a )c] 
[(R a + L a )< 



i a 



bSt dpb 

1 d 2 
f 2b 2 6t dpj 



[(Rb-L b )c] 



[(R b +L b )c]. 

(31) 



Using St = l/(a + b) , we identify the drift velocities for 
the two subgraph densities as 



a 

a + b 



Vh 



(Ra — L a ) — pb — p a , 
(Rb — Lb) = p a — Pb, 



(32) 



that again illustrates the early-time bias in the VM dy- 
namics on the complete bipartite graph. This two-time- 
scale dynamics is illustrated in Fig. 01 where a bipar- 
tite network of size a = b = 10 5 is initialized with 
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(p a (Q))Pb(0)) = (1,0). The dotted curve shows the con- 
vective transient that lasts for 3 time steps before the 
subsequent diffusive approach to the final consensus af- 
ter 9999.0 steps (solid curve). We define the end of the 
transient when the trajectory first approaches to within 
the extent of diffusive fluctuations about the diagonal. 

To determine the exit probability, we use the fact that 
lo = (p a + pb)/2 is conserved [see Eq. (fT5j) ]. so that 



1 



£l(Pa,Pb) = w(p a ,Pb) = TT^a + Pb)- 



(33) 



Strikingly, when one subgraph contains only Os and the 
other only Is, the probabilities of ending with all Is is 
1/2, independent of the subgraph size. To determine the 
average time T(p a , pb) to reach consensus — either all Is 
or all Os — as a function of p a and pb, we follow the same 
steps [lH as in Eqs. (|22|) and (|23|) to obtain the backward 
Kolmogorov equation 



NSt = (p. 



Pb) 



dT 
dp a 



dT 
dpb 



:(/Oa+Pb-2p a p b ) 



1 d 2 T 
a dpi 



\d 2 T 
bdpj 



(34) 



The first term on the right accounts for the bias that 
drives the system to equal subgraph densities while the 
second term accounts for the diffusion that ultimately 
drives the system to consensus. 

Since the subgraph densities p a and pb asymptotically 
approach each other and lo = (p a + pb)/2, we replace p a 
and pb by lo and also g g by \ in Eq. (fM)) to yield 



d 2 T 



(35) 



For the boundary conditions T(0) = T(l) = 0, the solu- 
tion is [compare with Eq. (|25p]. 



T(u) 



4ab 
a + b 



(1 - uj)ln 



1 



h lo In — 

LO LO 



(36) 



The consensus time has a similar form to that of the 
complete graph, but with an effective population size 
N e g = Aab/(a + b). If both components of a bipartite 
graph have a similar size, a, b rj N/2, then N e s ~ N. 
However, if the two components are disparate, e.g., for 
a - 0(1) and b w N then T - 0{l) ! Thus one highly- 
connected node strongly facilitates reaching consensus. 



B. The Two-Clique Graph 




0.2 0.4 0.6 0.: 

Pa 

1 W^-----^^** 




0.2 0.4 0.6 0.: 

Pa 




FIG. 4: (Color online) Subgraph densities p\>(t) versus p a (£) 
for single realizations of the voter model on (top to bottom): 
(a) a bipartite graph of 2 x 10 5 nodes; (b) a 2-clique graph of 
2 x 10 4 nodes, with (7 = 1 (upper trajectory) and C = 100 
(lower trajectory); (c) a Molloy-Reed graph of 2 x 10 nodes 
with degree distribution ~ k~ 2 5 . For the last example, 
the densities pa and pio are shown. 



Another instructive example that reveals the two-time- 
scale route to consensus is the two-clique graph (Fig. [5]). 
This graph consists of two separate complete graphs of N 
nodes, a and b, with each node in one clique also possess- 
ing C random cross links to the nodes in the other clique. 
As we now show, the extent of cross-linking strongly af- 
fects how the system reaches consensus. For the two-clique network, the respective probabilities 
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FIG. 5: (Color online) The two-clique graph with N = 12 and 
C= 1/2. 



that the clique density p a increases or decreases are: 



1 



N 



C 



N + C 



N + C 



Ph 



N , C , 



N + C 



N + C 



(37) 



with analogous expressions for the rates on clique b. For 
R a , the leading factor i is the probability of choosing 
clique a and (1 — p a ) is the probability of choosing a 
in clique a. The density of Is in clique a increases if 
we choose a in clique a (factor N ^ c p a ) or in clique b 

(factor N Z^ C Pb)- Similar explanations apply for L a and 
for the operators on clique b. 

The drift velocity for the density p a is then 



Spa. 
St 



(R a — L a ) 



C 



N + C 



( P b - Pa) • (38) 



When C is 0(1), the drift term is 0(1/N) and is therefore 
comparable to the diffusive term, 



D » = 2 — {R * + La)0C ^/N-N 



N 



C 



(39) 



which remains 0(1/ N) for all densities. The clique den- 
sities /5 a and pb thus evolve diffusively and indepen- 
dently until consensus which is reached. This behavior 
is illustrated in the upper dashed (green online) trajec- 
tory of Fig. Q|b) for a single realization of a 2-clique 
graph with a = b = 10 4 nodes in each clique and 
with a single crosslink per node for the initial condition 

(fla(0),pb(0)) = (0,l). 

Conversely, for C > 0(1) the drift is dominant and 
the clique densities are driven to a common value that 
subsequently evolves diffusively along the p & = pb diag- 
onal until consensus is reached as illustrated in the lower 
trajectory in Fig. [U[b). Shown is the trajectory of a sin- 
gle realization when there are 100 crosslinks per node for 
the initial condition (p a (0), pt>(0)) = (1,0). The dotted 
(blue online) line shows the initial transient of 100 time 
steps and the solid (red online) line shows the subsequent 
diffusive approach to consensus at p a = Pb = 1 in 6437 
time steps. 



C. Heterogeneous-Degree Networks 



transition probabilities for increasing and decreasing the 
density of voters of type 1 on nodes of fixed degree k are: 



R k [{p k }]=P[p k -^p+} = - ^'i- A xy $(?/, x) 



1 v^' 1 

L k [{p k }}=P[p k -» p k \ = — — A xy $(x, y) 

x,y 



(40) 



where p^ = p k ± N^ 1 , and the prime on the sums again 
denote the restriction to nodes x with fixed degree k. In 
Eq. (|40|) the densities associated with nodes of degrees 
k' ^ k are unaltered. 

We now make the approximation that the degrees 
of neighboring nodes are uncorrelated, such as in the 
Molloy-Reed (MR) network [13] ■ Then we may replace 
the elements of the adjacency matrix elements by their 
expected values to give 



.4 



xy 



(A X y) — 



kxky 



(41) 



This relation expresses the fact that in the absence of 
degree correlations, the probability that nodes x, y are 
connected is proportional to k x k y , and the proportion- 
ality constant in (|4"Tj) is determined by using the fact 
that the average node degree is just p,± = ^2 X A xy 
Substituting the mean-field assumption (|41[) for A xy in 
Eq. gO]), and using Eqs. (J), TO; and CH, tnc transi- 
tion probabilities simplify to 

Rfc = njtu;(l - pk), Lfe = n k p k (l - u). (42) 

Similar to Eq. (|22p , the recursion formula for the mean 
consensus time, starting with initial densities {p k }, is 

T({ Pk }) = 6t + J2n- k [{p k }}T(p+) + J2lk[{Pk}}T( Pk -) 

k k 

+ [l-J2(K k [{Pk}] + L k [{p k }})]T({p k }). 

k 

(43) 



Expanding (143)) to second order in N k gives the back- 
ward Kolmogorov equation for the consensus time 



V- dT V- r, 9 T 

V dpk k dPk 



(44) 



with degree-dependent velocity and diffusion coefficients 



- 6pk fv> t ^ 
v k = -^( R fc - Lfe) = w - p k , 

D k ^^(K k +U) = UJ + Pk ~ 2UJPk 



(45) 



2St 



2Nn k 



We now study the VM on networks with arbitrary de- 
gree distributions. Following Eq. (|3|), the fundamental 



To simplify the backward equation (|44p for the consen- 
sus time, we start with the forward Kolmogorov equation 
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for the probability distribution 

5 k , t) 

B fe [pfe + Sk]c(pk + 5k, t) 

k 

[l-^(F fe [{p fc }]+B fc [{p fc }])]c({p fc },t), 

k 

(46) 



c{{pk},t+5t) = ^¥ k [pk - 5 k ]c(p k 



+ 



and expand to second order to obtain the Fokker-Planck 
equation 

d d 2 

d t c({ Pk },t) = -Y,g-VkC + Y,g^ D kC- (47) 
k k ^ k 

Now we compute the mean value of the average density 
Pk with respect to the above probability distribution to 
obtain the time dependence 



d{pk) 
dt 



(to) - {pk}- 



whose solution is, using the conservation of to, 
(pk(t))=Lu(0)-[Lo(0)-p k (0)]e- t . 



(48) 



(49) 



Thus after a time scale that is of the order of one, all the 
Pk approach the common value of the conserved quantity 
lj and the drift velocity in (I45[) vanishes. This dynamics 
is illustrated in Fig. |3{c) where the trajectory of a single 
VM realization is shown for a Molloy-Reed network of 
2 x 10 5 nodes, with degree distribution nk ~ fc~ 2 ' 5 . The 
initial state is (pk>in (0), Pk<m (0)) = (0,1). Shown are 
Pe(t) (degree less than p\ = 8) and pio{t) (degree greater 
than pi) versus to. The initial transient, which lasts 1.78 
time steps, is shown dotted, while consensus occurs after 
1742 time steps. 

As a result of this rapid approach to the locus pk = <-o 
for each k, we may drop the drift term in Eq. (|44[) and 
also convert to derivatives with respect to to by using 



dT 
dpk 

to reduce ([44]) to 

-r 

k 



dT duj 
doj dpk 



knk dT 
pi dio ' 



Mi / 



dto 2 



We now define the effective population size by 



N 



Mi 

P2 



(50) 



(51) 



and comparing Eq. ([5T?]) with (|24[) . the consensus time is 



T N (w) = N eS 



(1 - w)ln- 



1 



CO 



a; In — 

10 



(52) 



To compute N c g for a network of N nodes with a 
power-law degree distribution, ~ k~ v , we first deter- 
mine the maximum degree in the network. We obtain this 
quantity from the extremal criterion jla |. J k k~ v dk — 

■if. This condition gives fc max ~ N 1 ^ v ~ 1 \ which holds 
for all v > 2, while for v < 2, the maximum degree re- 
mains of the order of N. With this upper cutoff, the m th 
moment of the degree distribution is 



Mm 



J k m n k dk 

'jY(m+l-^)/(^l) 

In TV 
1 0(1) 



m + l-v > 0, 1/ > 2 
m+l-v > 0, i/ < 2 
m+l—v = U 
m + l — v < 0. 



Using these results, the asymptotic behaviors of N e s and 
Tjv are 



(N 

N/\nN 

[0(1) 



v > 3, 
v = 3, 
2 < < 3, 
i/ = 2, 
z/ < 2. 



(54) 
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FIG. 6: (Color online) Consensus time Tjv versus N on the 
Molloy-Reed network with degree distribution nt = k~ v for 
v = 2.1 (+), 2.3 (x), 2.5 (*), 2.7 (o) and 2.9 (•). Each data 
point is based on 100 realizations of the graph and 10 realiza- 
tions of the voter model on each graph. The lines represent 
the theoretical prediction of Eq. (|54[1 . The inset shows the 
same data plotted in the scaled form ^Tn/pI versus N. 

The main feature of Eq. (|54|) is that consensus is 
achieved quickly for the VM, that is, Tjy -C N for all 
v < 3 (l6| . This consensus time is also much faster than 
the corresponding behavior on regular lattices in d spatial 
dimensions [III]: T N oc N 2 for d = 1 , T N oc N In N for 
d = 2, and Tjv oc iV for d > 2. We tested our prediction 
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(|54| by simulations of the voter model on both the MR 
network (Fig. [5]) and also a network that is grown by the 
the redirection algorithm [l9|]. This method generates a 
network with shifted linear attachment rate; that is, the 
probability of attaching to a node of degree k is given by 
k + A. This growth rule leads to a power-law degree dis- 
tribution network with tunable exponent v = 3 + A. For 
a network that is generated by the redirection algorithm, 
the N dependence of Tjv is essentially the same as that 
for the MR network (see Fig. 3 of Ref. [Hj]), and both 
these numerical results are in good agreement with our 
theory. 

One important feature of the network that is built by 
the redirection algorithm is that degrees of neighboring 
nodes are correlated In spite of this correlation, 

the actual values of the consensus times for the MR and 
the redirection networks are numerically within 15% of 
each other. Thus evidently degree correlations have a 
secondary role in determining the consensus time; the 
broadness of the degree distribution is much more im- 
portant. 



V. INVASION PROCESS 

We now study the role of degree heterogeneity on IP 
dynamics, a question that was first studied by Castellano 
6] . In analogy with our approach for the voter model, we 
replace the adjacency matrix elements by their average 
values, as in Eq. (|4Tj) , so that the transition probability 
of Eq. (jU becomes 



From this expression, and following exactly the same 
steps that led to Eq. (|12"|) , the transition probabilities 
for nodes of a fixed given degree are: 

k k 
Rfe = n k — p(l-pk) L fe = n fe — p k (l-p), (56) 
Mi Mi 

where again p^ = p ± N^ 1 . The corresponding fc- 
dependent drift velocity and diffusion coefficient are then: 



or mi 



(Spk 



2St 



2Nn k p-\ 



(57) 



From the drift velocity, the time dependence of the aver- 
age density p k is determined from 



d (Pk) k 

—rr = —{(P-Pk))- 
at mi 



(58) 



which shows that the densities p k again approach the 
common value p after a time scale of 0(1). Once this 
concurrence happens, we may asymptotically replace p, 



and concomitantly all the p k , by ui-i in the expressions 
([57)) for v k and D k . 

With this simplification, we now follow exactly the 
same steps that led from Eq. (|43|) to Eq. (|52|) in the pre- 
vious section to write backward Kolmogorov equation for 
the consensus time: 



LU-l (1 - LU-l) 



d 2 T 
dui-i 



(59) 



By comparing with Eqs. (|52|) and (|59| . we deduce the 
average consensus time 



T/v(w_i) = N cS 



(1 -w_i)ln 



1 , 1 

h U)-i 111 

1 — LJ-X U)-l 



(60) 

with iVoff = Npip-i. For graphs with power-law degree 
distributions we use the moments written in Eq. (|53|) to 
obtain 



T N tx N eS 



N v>2, 
NlnN v = 2, 
N 3 -» v<2. 



(61) 



Thus for all reasonable graphs, those with v > 2, the 
consensus time for the IP is strictly linear in N, in sharp 
distinction to much faster consensus for the VM. 



VI. DYNAMICS WITH SELECTION 

We now study the VM and IP when the two states 
and 1 have different fitnesses. We define resident Os 
as having fitness / = 1, and mutant Is having fitness 
f = r, which can, in principle, be either smaller or larger 
than 1. Our goal is to determine the fixation probability, 
namely, the probability that a single fitter mutant (i.e., 
we consider only the case r > 1) overspreads a population 
under biased dynamics [2(| ■ Such a phenomenon provides 
a natural description for epidemic propagation (2ll . [221 . 
l23l| . th e emergence of fads 0, HH, [2y| , social cooperation 
[8ll27ll28fl. or the invasion of an ecological niche by a new 
species [1, H, [l(| • 

We define the update steps in the biased VM as: 

(i) pick a voter with probability proportional to its in- 
verse fitness (f~ 1 /N{f- 1 )); 

(ii) the voter adopts the state of a random neighbor. 

Here (/ _1 ) is the mean inverse fitness of the population. 
Thus a weaker voter is more likely to be picked and to 
be influenced by a neighbor. We can equivalently view 
the fitness as the inverse of a death rate. Similarly, the 
evolution steps in the biased IP are: 

(i) pick an invader with probability proportional to its 
fitness (//JV(/»; 

(ii) the invader exports its state to a random neighbor. 
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We denote by (/) the mean fitness of the population that 
we can equivalently view as the number of offspring pro- 
duced by an individual. Thus a fitter mutant is more 
likely to spread its progeny. For both the biased VM 
and biased IP, we focus on the weak selection limit, with 
r = 1 + s and s<l. 

Thus far, we have studied the VM and IP in discrete 
time, where the time is incremented by 1/N after each 
update event. For the biased models, it turns out to be 
simpler to treat continuous-time versions. Thus for the 
biased VM with continuous dynamics, each individual 
adopts the state of a random neighbor at a rate propor- 
tional to its inverse fitness. Similarly, in the biased IP, 
each individual exports its progeny to a random neigh- 
bor at a rate proportional to its fitness. This rate-based 
update rule is equivalent to a discrete-time update with 
the time increment chosen from an exponential distribu- 
tion, with mean value 1/N(f~ 1 ) for the voter model and 
1/N(f) for the invasion process. The discrete and con- 
tinuous models arc then equivalent except for this overall 
factor in the time scale. 



A. Biased Voter Model 

For the biased VM on a general network, the proba- 
bility of changing the state at node x is, in close analogy 
with the unbiased case [compare with Eq. 



A, 



Nkxif- 1 ) 



$(y,x) + -$(x,y) 



, (62) 



Following the same approach as in Sec lIVC] the density 
Pk of Is at nodes of degree k increases by 5p k = 1/N k 
with probability Rfe(7?) and decreases by 1/N k with prob- 
ability Lfc(ry) in a single update, where 



Rk{v) = Nk(f-i) T,' A *y®(y> x ) 



3'!/ 



Mv) 



(63) 



are the respective transition probabilities for (0 — » 1) and 
(1 — » 0) in the weak selection limit s <C 1. The primes 
on the sums again denote the restriction to nodes x of 
degree k. We seek the fixation probability £ i to the state 
consisting of all Is as a function of the initial densities 
of 1. This probability obeys the backward Kolmogorov 
equation B£i = [H|,[l4j], subject to the boundary con- 
ditions £\(p = 0) = and £\(p = 1) = 1. In the diffusion 
approximation, the generator B of this equation may be 
expressed as a sum of the changes in pk over all fc, 



B 4s 



^(R fc -L fc )— + — 



(R fc + L fe )— j 



with St = l/Nif- 1 ). 



if 

dpi 
(64) 



On degree-regular graphs, the sums in Eq. (|63| include 
all nodes and hence count the total number of active 
links. For uncorrelated node degrees, the fraction of ac- 
tive links reduces to p(l — p), and the generator for 
s<l becomes 



B « p(l - p) 



d_ 

dp 



Ndp 2 



(65) 



The drift and diffusion terms differ by a factor O(sN), 
so that selection dominates when the population size N 
is larger than 0(1/ s), while diffusion, or random genetic 
drift, dominates otherwise. Since the probability of in- 
creasing the density of Is at each update is r times larger 
than the probability of decreasing this density, the fixa- 
tion process is the same as the absorption of a uniformly 
biased random walk in a finite interval. The fixation 
probability is thus given by the well-known exact formula 

m 



1 _ r -Np 



1 



-N 



(66) 



In the small-selection limit s <C 1, this result ap- 
proaches the solution of the backward Kolmogorov equa- 
tion B£i =0, with generator B given in ([B^]) . 



Sx(p) = 



1 



-sNp 



1-e 



-sN 



F(sN,p). 



(67) 



where for later usage we introduce the notation T for the 
fixation probability. 

Now we return to degree-heterogeneous networks. 
Here the conserved quantity for unbiased dynamics is the 
average degree- weighted density u> [Eq. {H}]. This con- 
servation law suggests that we study the evolution of (u) 
when r ^ 1. When node x is updated, u> changes by 

lu(t] x ) - cj(r?) = k x (l - 2r ] (x))/p. 1 N, (68) 

where r) x again denotes the state of the system after the 
update. Thus (cu) evolves as 



piN 



Y^A xy $(x,y), 



(69) 



x.y 



for sCl. Now we use the mean-field assumption ((41 
for A xy to reduce Eq. (fBU)) to 



R*(»7) = "^piy n k uj(l - p k ), 



1 - s 



(70) 



,j^T\ n kPk(l - w), 



and Eq. (J69J) to 



at 



8(U)(1-(U)), 



(71) 
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whose solution is 

<^)) = {l-[l-c(0)- 1 ]e- s *}- 1 . (72) 
Then the time evolution of (pk) becomes 
d{pk) ^(Rfe-Lfc) 



dt 



st 



(u)-{p k )+a(l-{w)){p k ). (73) 
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FIG. 7: (Color online) Illustration of the two-time-scale dy- 
namics. Moments of the 1 density in the biased VM and 
the biased IP on a Molloy-Reed network of 10 4 nodes with a 
power-law degree distribution rife ~ k~ v and v — 2.5. Nodes 
with degree larger than the mean degree are initialized to 1 
while all other nodes are 0. Here s = 8.13 x 10~ 4 for the VM 



ands = 6.6 x 10"^ 
1/iVeff, with 7V cff 
for the IP. 



for the IP. These s values were chosen to be 
= Npi/fj.2 for the VM and N e « = iV/xi/Lt-i 



To solve this equation we combine it with Eq. (fTTj) to 
yield 



d(uj - p k ) 



Of 



-(1 - s)(uj - p k )(l -u), 



with solution 



(p k (t)) = ( W W}- e -*[ W (0)- Pfc (0)]{ W (0) + [l- W (0)]e- s *}. 

(74) 

For small selective advantage (s <C 1), this equation in- 
volves two distinct time scales. In a time of the order of 
1, all the pk become equal to u>, the conserved quantity 
for the unbiased VM. Subsequently, the evolution of u> 
itself occurs on a longer time scale of order s _1 3> 1 and 
is driven by the bias (Fig. [7]). 

We now determine the fixation probability by replacing 
the pk by uj in the transition probabilities R and L in 
Eqs. ((701), and the derivative by ^-7^. Then the 
generator 



becomes, in the limit s«l, 



B 



w(l 



/<2 



8 



2 1 







(75) 



which is the same as the generator for degree-regular 
graphs in Eq. (|65p with the population size N replaced by 



N c ff — Np\jpi. Consequently, the solution to B£i = 
is simply f 1(0; ) = J-(sN c s , lu), with J 7 defined in Eq. (|67[) . 
that is in excellent agreement with our simulation data 
(Fig.©. 




FIG. 8: (Color online) Scaling plot of fixation probabilities 
for VM and IP dynamics for a Molloy-Reed graph with degree 
distribution ~ and v = 2.5, with N — 10 3 and fii — 8. 
The empty symbols correspond to IP dynamics with s = 0.004 
(□), s = 0.008 (o) and s = 0.016 (A); the filled symbols 
correspond to VM dynamics with s = 0.01, (■), s = 0.02 
(•) and s — 0.08 (a). The smooth curve is the prediction of 
Eq. (HI). 

When a single mutant is initially located at a node 
of degree k, then ui — k/Npi. Substituting this value 
into ^(sNcHjLlj), we obtain the general result that the 
fixation probability starting with a single mutant on a 
node of degree k is proportional to k; £\ cx k for all s <C 1 
(Fig. [5]). More precisely, £\ = ^(sNcg,^ — k/Npi) has 
two distinct limiting behaviors: 



1 



£1 



spi 



k s < l/N eS ; 



(76) 



k 1/JVeff < S < 1. 

M2 



In the s 1/N c g limit, we recover the results given in 
Eq. (fTS]) for unbiased evolution. Notice also that the 
relative influence of selection and random genetic drift is 
determined by the variable combination siV G ff. Because 
iV e ff can be much less than N, diffusion can be important 
for much larger populations compared to degree-regular 
graphs. 



B. Biased Invasion Process 

In the complementary biased invasion process, each 
individual reproduces at a rate proportional to its fitness. 
Hence the transition probability is 



A X y 



Nk y (f) 



[r$(y,x)+$(x,y)} . (77) 
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the node degree, E\ ~ 1/k (fig. [9] The limiting behaviors 
of the fixation probability are 

£1 = { (8i) 

s 1 

— 1/JV < s < 1. 

In the s <C 1/iV limit we again recover the results in the 
neutral case (jT5j) . 

VII. DISCUSSION 



FIG. 9: (Color online) Fixation probability of a single mutant 
initially at a node of degree A: on a Molloy-Reed network with 
n k ~ k~ v and v = 2.5, with TV = 10 3 and fii = 8. The 
empty symbols correspond to IP dynamics with s — 0.004 
(□), s = 0.008 (o) and s = 0.016 (A); the filled symbols 
correspond to VM dynamics with s = 0.01, (■), s — 0.02 (•) 
and s — 0.08 (a). The solid lines, with slopes +1 and —1, 
correspond to the second of Eqs. (|76|l and (|81|) . 



For degree-uncorrelated graphs and weak selection s -C 1 , 
the transition probabilities for increasing and decreasing 
Pfe in a single update become 

(f + s)k 
tik(V) = rp- n k p(l - pk), 



k 

L k(v) = — 77T n k p k (l - p). 
Pi\J) 



(78) 



Following the same steps that led to Eq. ([75]) and using 
8t = 1/N(f) for continuous-time dynamics, p k evolves as 



dpk 
dt 



Mi 



[p - Pk + sp k (l - p)} , 



(79) 



which implies that all the p k rapidly become equal, and, 
as a consequence, all moments uj m also become equal. 
For the unbiased IP, the conserved quantity is ct>_i; thus 
for weak selection, u>—i is now the most slowly changing 
quantity, as illustrated in Fig. [7] Hence we replace all p k 
by W—x, and transform all derivatives with respect to p k 
to derivatives with respect to in the generator l|64p 
to obtain 



B 



PiP-i 







du>- 



i d 2 

Ndto 2 ~ 



(80) 



from which, in close analogy with our previous anal- 
ysis of the VM, the fixation probability is £i(w_i) = 
F(sN, u>-i), with T again defined in Eq. (|67|) . 

As a basic corollary, consider the fixation probability 
when the initial state consists of a single mutant at a node 
of degree k. Substituting W—x — l/fciV/x_x into J-, we 
obtain the general results that in the small selection limit 
s«l the fixation probability is inversely proportional to 



We developed a unified framework to investigate the 
dynamics of the voter model (VM), the invasion process 
(IP), as well as link dynamics (LD) on complex networks. 
When the relative fitness of the two states, and 1, are 
the same, the variable q defined by 



VM, 

on regular graphs, 
IP, 



(82) 



is conserved by the dynamics. Because of this conserva- 
tion law, the exit probability £ 1 to 1 consensus is simply 
equal to this conserved variable. This simple result has 
far-reaching consequences on complex networks, where 
high-degree nodes can have a disproportionately strong 
influence. 

On complex networks, the evolution of the VM and 
the IP is governed by two disparate time scales. Ini- 
tially there is a quick approach to a homogeneous state 
in which the density of Is on nodes of degree k becomes 
independent of k for any initial condition. Diffusive fluc- 
tuations then drive the system to a final consensus state 
on a much slower time scale. For a network of N nodes 
with a power-law degree distribution, n k ~ k~ v ', the VM 
consensus time scales slower than linearly with iV 
for all v < 3. This result seems to hold independent of 
the extent of degree correlations. Thus when the degree 
distribution is sufficiently broad, high degree nodes facil- 
itate the approach to consensus. In contrast, in the IP, 
T/v scales linearly with N for power-law degree networks 
for all v > 2. 

We also studied the VM and IP with the additional 
feature of selection, in which individuals in state 1 are 
fitter than those in the state. The populations in these 
two models can be characterized by an effective size 



N, 



cff 



P2 



VM, 



(83) 



_Npip_i ~ N IP (& regular graphs). 



The equivalence N c g ~ N applies for the physically ac- 
cessible case where the mean degree of a graph is fi- 
nite. From the fixation probabilities given in the previ- 
ous section, the fundamental parameter for both models 
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is a = siVeg, which is just the Peclet number in the lan- 
guage of biased diffusion. The interesting limit is s — > 
and N — > oo, such that a is constant. Now the general 
form of the backward Kolmogorov generator is 



B 



JVcff 



9 a 2 
t — i 

dq dq 2 



(84) 



An important characteristic of the dynamics is the 
probability that a single 1 mutant overspreads an other- 
wise homogeneous population of Os. This fixation prob- 
ability satisfies B£x(g) = 0, with solution 



Si(q) = 



1 - e- aq 
1 - e~ a 



(85) 



that depends on the scaling combination sNeS and is in- 
dependent of s and N e g separately. Another important 
feature of fixation is its dependence on the degree of the 
node at which a mutant first appears. For the VM, the 
probability for a mutant on a node of degree k to fixate 
is proportional to k [Eq. Q760]. while in the complemen- 
tary IP, the fixation probability is proportional to 1/fc 
[Eq. (|8Tj) ]. The origin of this behavior is simple to un- 
derstand. In the VM, a well-connected mutant is more 
likely to be asked its opinion before the mutant queries 
one of its neighbors. In the IP, a mutant on a high-degree 
node is more likely to be invaded by a neighbor before 
the mutant itself can invade. Thus network heterogeneity 
leads to effective evolutionary heterogeneity. 

The above results also allow us to understand the fix- 
ation probability for the case where a mutant sponta- 
neously appears at random in the network. For the VM 
in the selection-dominated regime (sN c g 3> 1), averag- 
ing Eq. ([76]) over all nodes gives a fixation probability 



that is smaller by a factor fi 2 / ' n% < 1 than that on reg- 
ular graphs. Thus a heterogeneous graph is an inhos- 
pitable environment for a mutant with VM dynamics. 
The source of this inhospitability is that a spontaneously- 
appearing mutant is likely to be on a low-degree node, 
and its state is then quickly erased by interactions with 
higher-degree nodes. Conversely, for the IP, averaging 
Eq. (fSTj) over all nodes gives a fixation probability that 
is independent of the node degree. 

Finally, with our general formalism, we can determine 
the average consensus time from the solution of the back- 
ward Kolmogorov equation BT(g) = — 1. For the un- 
biased VM and the unbiased IP, this solution has the 
generic form: 



T N (q) = N cS 



{1-q) In 



1-9 



q In • 



(86) 



where q is the conserved quantity given in Eq. (|82p . In 
the biased case the solution of Eq. (|84l) can be written 
as an integral but it is quite cumbersome [1 21 ] . From the 
form of the generator in Eq. ([84]) , however, it is clear 
that T^(q) = iV e ff r(<7, a), and r(q,a) is a function that 
is non singular at a = 0. Hence, the dependence on 
effective population size N e s is not affected by the bias 
when a is kept constant. 
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